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Abstract 

The eigenvalues of the transfer matrix in a six-vertex model (with periodic bound- 
ary conditions) can be written in terms of n constants vi,...,v n , the zeros 
of the function Q(v) . A peculiar class of eigenvalues are those in which two 
of the constants v\,V2 are equal to A, — A , with A = —cosh A and A re- 
lated to the Boltzmann weights of the six-vertex model by the usual combina- 
tion A = (a 2 + b 2 — c 2 )/2ab. The eigenvectors associated to these eigenvalues 
are Bethe states (although they seem not). We count the number of such states 
(eigenvectors) for n = 2, 3, 4, 5 when N , the columns in a row of a square lat- 
tice, is arbitrary. The number obtained is independent of the value of A , but 
depends on N . We give the explicit expression of the eigenvalues in terms of 
a,b,c (when possible) or in terms of the roots of a certain reciprocal polyno- 
mial, being very simple to reproduce numerically these special eigenvalues for ar- 
bitrary N in the blocks n considered. For real a, b, c such eigenvalues are real. 

PACS numbers: 05.50. +q 75.10.Hk 
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1 The problem 

Some time ago the author of this note read in the paper Completeness of the Bethe Ansatz for 
the Six and Eight-Vertex Models by R.J Baxter [1, Sect. 4] the following sentence concerning 
certain proper states of the transfer matrix in the six-vertex model at zero-field: 
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The other problem that we encountered first occurs for N — 4 and n = 2 , then for even N 
and 2 < n < N — 2 . It is referred to by Bethe himself and has been considered by others since 2 . 
For some eigenvalues with momentum ±1 , 3 i.e. k\ + ■ ■ ■ + k n = or tt , we found that 

n 

Q(v) = Y[smh[(v- Vj )/2] 

3=1 

had a pair of zeros vi 7 v 2 such that v\ = A, v 2 = —A . 

The lines continued later as follows: 

For N = 4 there was just one such eigenvalue A , in the n = 2 central block. For N = 6 there was 
one in the n = 2 block, two in the n = 3 block, and one in the n = 4 block. For N — 8 there were 
1,2,5,2,1 in the n = 2,3,4,5,6 blocks, respectively. This suggests (tentatively) that the Catalan 
numbers may count such eigenvalues. 4 The momenta were —1 except for a single eigenvalue with 
momentum +1 in each block with 3 < n < N — 3 . 

If the author had understood properly the eigenvectors associated to such eigenvalues and 
how to obtain them from Bethe ansatz, probably would have not detained so long when reading 
these sentences. But that was not the case: we were calculating at that time the free-energy 
per site of a vertex model whose ground state was a state of this type, and the value of the free- 
energy that we were deriving was once and again the incorrect one. We decided in consequence 
to put aside the free-energy problem for a time and study instead these states in the six-vertex 
model. We ignore the correct name that we shall use for them. In the literature they have 
received the name of singular Bethe states or singularities of the Bethe solutions [3, 4, 6], and 
also non-Bethe eigenvectors [7]. We might even remember some references in which they are 
alluded as improper states. Since they need a name and no other states are considered in this 
paper we will refer to them as bound pairs merely. 

This note communicates some results of the study and answers the interrogation suggested 
in Baxter's paper: Are Catalan numbers counting the bound pair states of a square six-vertex 
model with periodic boundary conditions? 

2 The model 

The model to be considered is a six- vertex model in a square lattice [8, 9]. In this model to 
each site of the lattice is associated one of the six arrangements of arrows shown in figure 1, 
where each of these arrangements has an energy e\, . . . , Eq and a Boltzmann weight given by 

tuj = exp(-£j/k B T), j = l,...,6. 

The configurations of arrows satisfy the 'ice rule', because at each site of the lattice there are 
two arrows in and two arrows out. 

2 [2, after eq. (23)], [3], [4], [5] 

3 Baxter means e*( fcl H h fen ) = ±1 

4 Catalan numbers are 1, 2, 5, 14, 132, 429, . . . 
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1 2 3 4 5 6 

Figure 1. The six configurations allowed at a vertex. At each site of the lattice there are two arrows in and two arrows 
out. This is known as the 'ice-rule'. 

Suppose that the lattice has dimensions M x N , that is N sites horizontally and M verti- 
cally, with the imposition of periodic boundary conditions in both directions. The state of an 
arbitrary row of N vertical edges is then specified by the configuration of up and down arrows 
on the edge. Let er = (a±, . . . , <jn) denote the state ( aj = +1 for an up arrow at vertex j , 
Gj = —1 for a down arrow). If a is the state of a row and a' the state of the row bellow, 
the two adjacent states are coupled by the transfer matrix T aa / , whose entries are given by a 
trace of 2 x 2 matrices 



T au , = trace R aia ^ R a , 2a > 2 



R 



(2.1) 
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R 



++ 



LOl 

0^4 



R 



+- 



uj 5 




R. 



UJQ 





R. 



uj 3 

Ll>2 



A consequence of the 'ice rule' together with the horizontal periodicity of the lattice is that 
the number n of down (or up) arrows in a row is a conserved quantity from row to row, and 
T , a 2 N x 2 N matrix, breaks up into N + 1 diagonal blocks with one block for each value 



n 



0, 1, . . . , N . The dimension of block n is 



. The transfer matrix is used to calculate in 



statistical mechanics the partition function of the lattice Z = trace T M , and this has implied 
the diagozalization of matrix T . In the case of a zero electrical field (the case treated here) 
where 

a = 0J\ = Ld2, b = LO3 = UJ4, c = UJ5 = UJq, 



the eigenvalues A of the transfer matrix are known to be [9] 

(f>{\ - v) Q(v + 2A) + <P(X + v) Q(v - 2A) 



where functions 



a(v) = (-iy 

(v) , Q(v) are 



q(v) 



<f,( v ) = p N smh. N (v/2) 

Q(v) 



\ J sinh [(■? 
j'=i 



«i)/2], 



and p , X , v are defined so that 



(2.2) 
(2.3) 

(2.4) 
(2.5) 



a = p sinh - (A — v), 



b = p sinh - (A + v), 



c = p sinh A. 



(2.6) 



To write the eigenvalues (2.3) we have to locate v\,...,v n for these eigenvalues. There are 
many solutions, corresponding to the different eigenvalues. 
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3 Catalan numbers or not 

Do Catalan numbers count the bound pair states of a square six-vertex model with periodic 
boundary conditions? The answer is no. 

As was communicated in ref. [1], it is true that for N = 4 there is one bound pair in the n = 2 
central block, that for N = 6 there are 1,2,1 in the n = 2,3,4 blocks, respectively, and 
that for N = 8 there are 1,2,5,2,1 in the n = 2,3,4,5,6 blocks. However, if the counting 
started in the previous reference had continued it would have found that there are 1,2,6,10 
in n = 2, 3, 4, 5 5 for N = 10 , and 1,2,7,12 in n = 2,3,4,5 for N = 12 . In fact our 
calculations here show that for general even N the number is exactly 1, 2, N/2 + 1,N in the 
blocks n = 2, 3, 4, 5 . The numbers of states for n = 6 and beyond wont be studied in this 
paper. 



4 Bound pairs and Bethe Ansatz 

To obtain the eigenfunctions of the transfer matrix one can either diagonalize exactly the matrix 
(impossible when the size is not reasonable) or use the Bethe ansatz, the trial form that Bethe 
used for diagonalizing the quantum-mechanical Hamiltonial of the one-dimensional Heisenberg 
model [2]. The ansatz suggests that the eigenstate of T[y) , T(v)\ip) = A(v)\i/j) , can be written 
as \ip) = J2 Xl <...<x n f( x i> • • • • • • > where the coefficients f(x±, . . . ,x n ) are 

/(*!, ...,x n ) = J2 A Pl ,..., Pn e ik ^ ■ ■ ■ e ik ^. (4.1) 
p 

The numbers x±, . . . , x n indicate the positions of n down arrows on the lower vertical edges of 
a row of the lattice, and are ordered so that 1 < X\ < xi < . . . x n < N . We have experienced 
that the coefficients A pi) ... )Pn suitable to construct bound pairs 6 are given by 

A pi ,..., Pn = ep/C 11 s p . iPi , (4.2) 

l<i<j<n 

where ep = ±1 is the sign of the permutation {pi, . . . ,p n } of {1, . . . , n} and C is a non-zero 
constant to be fixed later in the most convenient manner (usually normalization). The vertex 
model defined by activities a, b, c so that 

A = a2 + b \~ c2 (4.3) 

2ab v ; 

enters in Sij , defined as 

Sij = 1 - 2 A e ik i + e i{ - ki+k i\ (4.4) 

To write the eigenstates is only necessary then to know the factors e lkl , . . . , e lkn that appear 
in (4.1) and (4.4). These factors are the solutions of the equations 

e Pl A P2 ,..., Pn , pi = A pii ... )Pn , (4-5) 



5 We omit to mention the blocks n = N/2 + 1 to N — 2 since the number of these states is the same as in 
the blocks 2, . . . , N/2 — 1 but in reverse order 

6 Because they give the same result that when the transfer matrix is directly diagonalized. When the size 
of the matrix allows it it is possible to carry many numerical experiments and they confirm this choice for the 
amplitudes 
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that impose the periodic boundary conditions on the problem making that f(x±,X2, • • • , x^) = 
f(x2,---,xpf,xi+N) what identifies the N+l and 1 vertices. To be (4.5) consistent equations 
among themselves, it is necessary that 

e iN(h + ■ ■ ■ + kn) = L ( 4 _g) 

Equations (4.1), (4.2), (4.4), (4.5) and (4.6), are sufficient equations to write bound pair eigen- 
functions, and when needed we will refer to them as "the Bethe ansatz equations for bound 
pairs". However, and this is not less important, it is also necessary a correct normalization 
of the eigenfunction. Without it, the state cannot be obtained. We have learned the correct 
normalization in ref. [1, Sect. 4], and show an example later for N = 6 and n = 3 . Equations 
(8.2), (8.3) and (8.2), (8.3), (8.6) are deduced taking into account such normalization. 

It is important to write, before finishing, the relation between k\ . . . ,k n and v\,...,v n in 
(2.5) (or better between e tk J and e v ' ) 

as mentioned in many papers. This relation permits to move from the eigenvalue (2.3) to the 
eigenvector (4.1) of the transfer matrix when we precise it. 



5 A change of variables 

Before describing any eigenvalue we make a useful change of variables concerning v and A in 
(2.6). The change is convenient for those (the author in this specific problem among them) 
who prefer to work with polynomials rather than with hyperbolic functions as in (2.5). Define 
the variables 

z = e~\ y = e- x (5.1) 
instead of v and A, then (2.5) is essentially 7 the polynomial in z and \jz given by 

1 n 

with Zj = e~ Vj , j = 1 . . . , n . To be correct we should have defined another symbol for (5.2), 
Q(z) for instance, however we will use the same letter with the understanding that Q(v) stands 
for (2.4) and Q(z) for (5.2). In terms of these variables and together with definitions (2.4) 
and (2.5), relation (2.3) becomes 

(2/p) N A(v) Q(z) = ±> [(z - yf Q(zy 2 ) + (1 - zyf Q(z/y 2 )] , (5.3) 

where multiplicative constant factors in Q cancel out of the calculations. To operate in a 
computer we prefer to work with this relation more than with (2.3). 

7 Essentially means up to multiplicative constants that do not depend on z (they may depend on y because 
y is regarded as a constant: after all y is fixed by the value that we choose for A , and viceversa). The constants 
are not relevant because do not change the value of A(v) , as commented after equation (5.3) 
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6 n=2 



This is the simplest case to study because the transfer matrix of N edges (with TV even) has 
only one bound pair state in this block for arbitrary A defined in (4.3). Since bound pairs are 
characterized by v± = A , v 2 = —A as mentioned in Sec. 1, function (5.2) factorizes as 

Q(z) = (zy-l)(z-y)/z, (6.1) 

the zeros of Q(z) being z\ = y and z 2 = 1/y. Introduced this function in (5.3) and noting 
that the r.h.s. is exactly divided by Q(z) in the l.h.s, the quotient affords the eigenvalue 8 

A(v) =a 2 b 2 (a N -* + b N ^)-c 2 (a N - 2 + b N - 2 ), N > 4, n = 2, (6.2) 

that is valid for generic N even. It can be checked numerically that (6.2) is always an eigenvalue 
of the transfer matrix for all values of a, b, c real or complex, 9 and since the block n = 2 is 
among the blocks of smallest dimensions, it can be done even for N not too small. The 
eigenvector associated to (6.2) was known to Bethe himself [2, also after eq. (23)] and is 
proportional to 

N 

|^} = ^(-l) i |M + l), (6.3) 
l=i 

after appropriate normalization. We do not reproduce here this eigenvector with the Bethe 
ansatz (the example that we reproduce is for N = 6 , n = 3 later), but want to comment 
about e tkl and e tk2 . The product of these two factors is for the eigenfunction (6.3) equal to 
— 1 , since from (4.1) derives the relation 

/(xi + l,x 2 + l) = e i ( kl + k ^f(x 1 ,x 2 ), with N+l = l, (6.4) 

which is simply a consequence of the translation invariance of the transfer matrix (2.1). But 
also v\ = X in (4.7) fixes e tkl = , what obliges to set 

e ih = _ e -ik 2 = Q) (65) 

as was done in [1]. This happens for all bound pairs that we have obtained no matter the 
values of N and n : it is simply a fact that for these states in this model 

e i(ki + k 2 ) = _ 1 (66) 

This condiction, together with the two identities in (6.5) mark how to work appropiately with 
bound pairs. 

8 A(v) in (5.3) is obtained in terms of z and y , of course. We have reexpressed the result in terms of a, b, c 
to write (6.2) 

9 a,b,c , the Boltzmann weights (2.2) of the vertex model, are real and positive, but when diagonalization of a 
matrix is considered in general, with no restriction to physical values only, they can also be negative or complex 
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7 n=3 



The trial function (5.2) is now of the form 

Q(z) = (zy-l)(z-y)(z- A)/z 3 ' 2 , (7.1) 

with A a constant (numerical or depending on y ) to be determined. Substituting (7.1) in (5.3), 
the r.h.s. of this equation is exactly divided by Q(z) in the l.h.s. if and only if ^4 = 0, —1, 1 
or A is the solution of a certain polynomial whose coefficients depend only on A . The root 
A = is not an admissible solution because (7.1) has not the required expansion (5.2); on 
the contrary, roots A = —1,1 yield admissible functions Q(z) because the associated A(i>) 
by (5.3) are always in the spectrum of the transfer matrix, as we have verified in numerous 
experiments. For example, the numbers 

A + = 2 a 3 b 3 - abc 2 (a 2 + ab + b 2 ) + c 4 (a 2 - ab + b 2 ), (7.2) 
A_ = 2 a 3 b 3 - abc 2 {a 2 - ab + b 2 ) - c 4 (a 2 + ab + b 2 ), (7.3) 

are eigenvalues of the N = 6 transfer matrix for arbitrary values of a,b,c. The first is for 
A = — 1 , the second for A = 1 . We present some of these numerical tests in Table 1. Regarding 
the situation in which A is the solution of a certain polynonial, when N = 6 such polynomial 
is 

A 4 + (8A 3 - 4A) ,4 3 + (20A 2 - 14) A 2 + (8A 3 - 4A) A + 1 = 0, (7.4) 

but it has to de discarded because none of the four roots of (7.4) is linked to an eigenvalue of 
the transfer matrix for arbitrary A (it can be checked also with Table 1). There are only two 
Q 's (that is, two bound pairs in the block) and two eigenvalues. 



N = 6, ra = 3 



z — exp (i 7r/5) 


z = exp (i 27r/5) 


z = 64 


y = exp (i 7r/3) 


y — exp (i 37r/4) 


V = 25 


p = i 


p = i 


p=l 


A = -1/2 


A = y/2/2 = 0.707... 


A = -313/25 = -12.52 






xl0 b 


0.779 508 


2.028 398 


24.899 897 


0.270 722 


0.421 472 


10.523 295 + 


0.168 519 


0.305 370 


8.909 579 


0.168 519 + 


0.074 622 2 


3.719 286 


-0.105 922 


0.065 240 7 


-3.724 777 


-0.229 696 


0.001298 21 + 


-8.840 959" 


-0.465 600 - 


-0.349 829~ 


-10.587 344 


-0.527034 


-0.448 262 


-24.913 779 


0.265 632 ±0.336 454 i 


0.851 314 ± 0.535 651 i 


12.449 161 ±21.567 683 i 


0.086 642 ±0.196 887 i 


0.395 856 ±0.293 548 i 


4.891 466 ± 7.822 856 i 


0.084 086 1 ± 0.159 490 i 


0.184 229 ±0.002 347 62 i 


4.750 077 ±8.912 538 i 


-0.084 223 2 ±0.124 642 i 


0.098 987 5 ±0.098 842 6 i 


-4.756 978 ± 8.869 031 i 


-0.131 171 ±0.269 132 i 


0.039 950 5±0.112 727i 


-4.887 688 ± 7.863 040 i 


-0.176 703 ±0.299 057 i 


0.003 395 86 ±0.451 325 i 


-12.457 140 ±21.571 004i 


a = sin(7r/15) = 0.207... 


a = sin(77r/40) = 0.522 . . . 


a = 39/80 = 0.4875 


b = sin(47r/15) = 0.743 . . . 


6 = sin(237r/40) = 0.972... 


b = -1599/80 = -19.9875 


c = V3/2 = 0.866 . . . 


c = y/2/2 


c = -312/25 = -12.48 



Table 1. In vertical are shown the 20 eigenvalues of the transfer matrix block iV = 6,n = 3 for different values of 
a,b,c . The eigenvalues are obtained by numerical diagonalization of the matrix in (2.1), and each result approximated 
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to the number arrayed in the table with the rule of 5 . In all the examples we have fixed z,y,p, and a,b,c are 
derived from them through (2.6). The values marked with + and — coincide, no matter the number of digits of 
accuracy demanded in the computation, with the theoretical values (7.2), (7.3) obtained in this paper solving (5.3). 
In the third column it is necessary to multiply by 10 6 to obtain the correct eigenvalue. Notice that when A = —1/2 
the bound pair (7.9) is degenerated and the transfer matrix has another linearly independent proper state with the 
same eigenvalue a 6 + b b . This degeneration happens for all values of a,b,c and not only for the particular value 
listed here. 

The situation is the same for arbitrary N even: there are only two bound pairs in the block 
and the generalization of (7.2) and (7.3) is 

n N-5,uN-5\ / n N-3, h N- 



A_ = a 3 b 3 (a N ' e + b N ~ 6 ) - abc 2 [ a N ' 4 + b N ~ 4 - ab - — u ) - c 4 / *' " " 



A + = a 6 b 6 (a n ~ b + b 1 "-*) - abc 1 a JV " 4 + 6 iV " 4 + ab - + c 

a+b J \ a + b 

\ 4 /a 

a — b J \ a 

that correspond to 

Q+(z) = (zy-l)(z-y)(z + l)/z^ 2 and Q~{z) = (zy - l)(z - y)(z - l)/z 3 / 2 , (7.5) 

respectively. The quotients written in A± above are fictitious 10 because the divisions can be 
performed exactly giving as result polynomials in a, b, c with no denominators. 

Each eigenvector of the transfer matrix has associated a given Q(z) , we now calculate as 

an example the eigenvector associated to Q + in (7.5) for N = 6 using Bethe ansatz. 11 For 
such state the product 

e i (fci + k 2 + h) = ^ / 7 _ 6 n 

that can be justified in several manners: one, if the eigenvalue is known, (7.2) in this case, it is 
enough to set 6 = 0, a = c in the eigenvalue. The coefficient of c N is precisely e*^ lH hfcn ^ 
[9]; or two, evaluating (— l) n Q(zy 2 )/Q(z) at the point z = 1/y [1]. This gives also such 
product. Since the third zero of the function Q + is at 23 = — 1 , relation (4.7) indicates that 
e ifc 3 _ _j ^ th&t substituted in (7.6) gives the product e l<ykl+k2 ^ = — 1 , something that seems 
to be shared by all bound pairs of the model as we remarked in (6.6). For our pair holds again 
(6.5) what makes that the factor S21 vanishes according to (4.4). To obtain the correct bound 
pair state the rule is 12 : calculate the Sij that do not vanish (in the present case there are 
five of them) with (4.4), keeping only the dominant term as e lkl goes to zero, and calculate 
S21 with (4.5). In this manner, instead of writting 'S21 = 0' in the formulae, S21 takes the 
expression that vanishes most rapidly as e lkl goes to zero. This expression is 

S21 = 2A (1 + 2 A) e i(N-l)h t (7.7) 



while 



s 12 = 2Ae s 13 = l + 2A, s 31 = 1, 

S23 = e~ ikl , s 32 = (1 + 2 A) e~ ikl . (7. 



10 I.e., introduced by the author to make the expressions compact 

11 We insist on the words Bethe ansatz because some authors refer to bound pair states as non-Bethe states, 
and they are Bethe states 

12 We have taken this rule from [1, Sect. 4] 
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Note that the amplitudes obtained with (4.2) after the substitution of (7.7) and (7.8) do 
satisfy exactly equations (4.5), as expected. Take now N = 6 . Inserting the values 
(4.2) into (4.1) we find that for example, /(1,2,3) = -2 A (1 + 2A) 2 /C and f (1,2,4) = 
2 A (1 + 2 A) e~^ 1 /C . In the case N = 6 two more components are necessary to write the 
eigenvector, namely 

/(1,2,5) = -2A (1 + 2A) e~ ikl /C, /(1,3,5) = -6A (1 + 2A) /C, 

since the remaining components are deduced from these four with the generalization of property 
(6.4) to the case n = 3. Clearly /(1,2,4), /(1,2, 5) are the elements that grow most rapidly 
as e lkl vanishes and the sensible choice here is to take C so that /(l, 2, 4) = 1 . The result is 
the right eigenvector associated to Q + in (7.5) 13 

\i/>) =|1,2,4) + |2,3,5) + |3,4,6) + |1,4,5) + |2,5,6) + |1,3,6) 

- |1, 2, 5} - |2, 3, 6) - |1, 3, 4} - |2, 4, 5) - |3, 5, 6) - |1, 4, 6), (7.9) 

which coincides with the vector found in [3, eq. (22)] using different methods. 
8 n=4 and n=5 

There is no problem in repeating the same steps as in n = 3 to deduce the number of bound 
pairs when n = 4 or n = 5 . In fact introducing 

Q(z) = ( Z y-l)( z -y)( z 2 + Az + B)/z 2 (8.1) 

into (5.3), it is possible to find constants A and B so that the function A(v) is an eigenvalue 
of the transfer matrix block n = 4 for arbitrary a, b, c activities. However, we follow a different 
method in this section with the intention of obtaining a better trial function Q not as general 
as in (8.1): we solve directly Bethe ansatz equations (4.5) instead 14 . The equations are already 
solved for e\ and e2 (for brevity we will use from now the notation e\ to denote the number 
e tkl , C2 to denote e lk2 , and so on), since we know that e\ = , e2 = — 1/ei , with the product 
eie2 equal to —1 as a characteristic of bound pairs. It remains to solve for e 3 , e 4 in the case 
n = 4 , and for e 3 , e 4 , e§ in the case of n = 5 . And when resolving the same care about s ij 
has to be taken that when the eigenfunction (7.9) was constructed in the previous section: S21 
that vanishes has to be evaluated with (4.5), taking then the expression that vanishes most 
rapidly as e\ goes to zero, and the remaining with (4.4). With these remarks taken into 
consideration the equations to solve are 

(8.2) 

N > 8, (8.3) 

1,5 For general N the eigenvector is \ip) = X^iG'i 2 + 1,2 + 3) — |2, 2 + 2, 2 + 3}) . The state that accompanies 
to Q~ is |V) = J2fLi(^ 1 )\\ 1 , 1 + 1,2 + 3) + |2,2 + 2,1 + 3)) . It has some similary with (6.3) but in the block 
n = 3 

14 Once e lk3 , . . . , e tkn are found solving Bethe equations, we use (4.7) to write Q given by (5.2) 



/ 1 - 2Ae 3 \ / l-2Ae3 + e 3 e4 \ 

V e 3 -2A J \l-2Ae 4 + e 3 eJ ' 
( \ - 2Ae 4 \ / l-2Ae4 + e 3 e 4 \ 

V e 4 - 2A J Vl-2Ae 3 + e 3 eJ ' 



9 



_ 



in the block n = 4 , and 

N _ x = / l-2Ae 3 \ / l-2Ae 3 + e 3 e 4 \ / l -2Ae 3 + e 3 e 5 \ 
63 " V e 3 -2A ^ Vl-2Ae 4 + e 3 e 4 y Vl-2Ae5 + e 3 e 5 y ' l ' j 

ef - 1 = f 1 " 2A64 V 1 ~ 2 f 4 + 2 f 4 + 6465 1 , iV>10 (8.5) 

V e 4 -2A J \l-2Ae 3 + e 3 eJ V 1 - 2Ae 5 + e 4 e 5 / ~ y 1 

_ / l-2Ae 5 \ / l-2Ae 5 + e 3 e 5 \ / 1 - 2Ae 5 + e 4 e 5 \ 

V e 5 - 2A J VI - 2Ae 3 + e 3 e 5 / \1 - 2Ae 4 + e 4 e 5 / ' 1 ' ' 

when n = 5 . Remember that A is given by (4.3) and N is an even number. 

Consider the equations relative to n = 5 for a moment. Notice that if (e 3 ,e 4 ,es) is a 
solution of equations (8.4)-(8.6) for given N and A 15 , also (e 4 ,e 3 ,es) , the interchange of e 3 
with e 4 , is a solution; and also it is (e 3 ,es,e 4 ) . Equations (8.4)-(8.6) do not distinghish a 
solution from any of its permutations. It is for this reason that two solutions are considered 
the same if coincide up to permutations. 

There is another relevant property of the equations: if (e 3 , e 4 , es) is a solution, -^J 
is also a solution for the same N and A . This feature brings considerable insight into the 
resolution of (8.4)-(8.6). For example, if e 3 is in the solution so does l/e 3 , as this property 
establishes, therefore l/e 3 is one of the numbers in (e 3 , e 4 , es) . If it is equal to its inverse, e 3 
is 1 or — 1 , but if not, the inverse of e 3 has to be say, e 4 , and thus e 3 e 4 = 1 . The argument 
is repeated with e 4 to conclude that e 4 is 1 or —1 or the inverse of e 3 . Finally, it is the 
turn of es , that can be only ±1 and not the inverse of any other number because there are 
no more left numbers to be paired with. In conclusion: (e 3 , e 4 , es) are (1, 1, 1) , (—1, —1, —1) 
or (e 3 ,e 4 ,±l) , with e 3 e 4 = 1 . There are no more possibilities for arbitrary A. Something 
similar happens when n = 4: the only solutions (e 3 ,e 4 ) of (8.2), (8.3) with A arbitrary are 
(1,-1) or the combinations (e 3 ,e 4 ) that satisfy e 3 e 4 = 1. Obviously this is so because the 
two properties explained above, permutation and inversion, hold for equations (8.2), (8.3) as 
well 16 . 

Lemma 8.1 (n = 4) The numbers e 3 ,e 4 given by equations (8.2), (8.3) subject to the condi- 
tion e 3 e 4 = 1 , are the roots of the quadratic polynomial 

x 2 - (r + l/r)x + 1 = 0, (8.7) 

where r is, in turn, the solution of the polynomial of degree N with coefficients fixed by A 
given by 

r N - 3 Ar^ 1 + 2 A 2 (r^ 2 + r 2 ) - 3 Ar + 1 = 0. (8.8) 



Proof Very simple. Just substitute directly e 3 = r , e 4 = 1/r in (8.2) and write the relation 
that results. Zero solutions r = are not wanted 17 . □ 

15 A fixed though arbitrary 

16 Observe that for all bound pairs obtained so far the product ei • • • e n = ±1 , something already mentioned 
in [1] and [3]. The momentum of these states, the sum of the k 's, is therefore or 1 (mod 2% ) 

17 We want 6364 = 1 with e-.j and finite numbers. Therefore none of them vanishes. We do not want more 
special objects like the pair eie2 = —1 with ei = 
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Surprisingly, the polynomial in (8.8) has the same coefficients when N = 8 , say, that when 
N = 100 , only that in this case the coefficients are distributed according to a degree 100 . 
Equality (8.8) belongs to the class of reciprocal equations [10] because the coefficient of r N is 
the same as the independent term, the coefficient of r N ~ l the same as the coefficient of r , 
and so on. If R is a root of a reciprocal equation, so it is its reciprocal 1/R. This cannot be 
a surprise, merely it is an expected consequence of the second property of the Bethe equations 
remarked a few paragraphs above. 

Lemma 8.2 (n = 5) The numbers e^^e^e^ given by equations (8.4)-(8.6) with the additional 
requirement e^e^ = 1 , e$ = — 1 , are the roots of the cubic polynomial 

(x + 1) (x 2 + (r + l/r)x + l) =0, (8.9) 

where r is the solution of (for simplicity we write the polynomial when N = 10 ) 

r w + (5 A + 2)r 9 +2 (2 A + l) 2 r 8 + 2 (2 A + 1)(A + l) 2 (r 7 + r 6 + r 5 + r 4 + r 3 ) 

+2(2A + l)V + (5A + 2)r + l = 0. (8.10) 

This is a reciprocal equation too. When N is arbitrary, the polynomial that generalizes (8.10) 
is a polynomial of degree N : r 10 , r 9 , r 8 above change into r N , r^" 1 , r N ~ 2 , respectively, 
and r 7 + • • • + r 3 into r N ~ 3 + • • • + r 3 . Nothing else changes. With these directions we avoid 
to write the generalization explicitly. 

When the requirement is 6364 = l,e§ = 1 , the solution (e$, 64,65) of equations (8.4)-(8.6) is given 
by (a; — 1) (x 2 + (r + 1/r) x + l) = , i.e., e 3 = — r, e4 = — 1/r, e$ = 1 , with r the roots of the 
polynomial obtained changing r by — r and A by —A in (8.10). The polynomial thus obtained 
is generalized to other N 's with the directions explained in the previous lines. 

Proof The substitution of 6364 = 1 and e$ = — 1 in (8.6) gives no information because the 
l.h.s. of (8.6) reduces to a number and the r.h.s. to the same number. However, substituted in 

(8.4) (or in (8.5)) is obtained a relation between the sum e^-\ = e 4 H = u and A . This 

e 3 e 4 

relation depends on N and, for example, when N = 10 is given by 

ii 5 -(5A + 2) u A + (8 A 2 + 8 A - 3) u 3 - (4 A 3 + 10 A 2 - 12 A - 6) u 2 

+ (4 A 3 - 14 A 2 - 16 A + l) u + 2 (2 A - 1) (A 2 + 3 A + l) = 0. (8.11) 

It is hard to see any recurrence in this equation but if u is decomposed into a number and its 
inverse, i.e., as u = — (r + 1/r) , r is a root of (8.10), which is a much simpler equation than 
the previous one. The numbers e^ = — r , e± = —1/r , e$ = — 1 , are therefore roots of (8.9) 
with r given by (8.10) if TV = 10 . □ 

Now we count states. Starting with n = 4 , we have the state characterized by 
(ei, e2, es, e^) = (ei, — 1/ei, 1, — 1) obtained before Lemma 8.1. For this state e^e-ze^ = 1, 



11 



and Q and A are given by 



Q(z) = (zy- \){z - y){z 2 - l)/z 2 , (8.12) 
A = a 4 b 4 (a N ~ 8 + b N ~ 8 ) - a 2 b 2 c 2 (a N ~ (i + b N ~ 6 - 2 a 2 b 2 " ' ~ ' 



a 2 -b 2 



N-6_ b N-6^ / a N-A_ b N 



- 3qW ( J +c T j' (8 - i3) 

as deduced from (4.7), (5.2) and the relation (5.3). As in A± obtained in Sect. 7, the quotients 
in (8.13) are artificial, and the divisions can be performed exactly giving for A an homo- 
geneous expression of order N in a, b, c with constant coefficients. Regarding the solution 
(ei, — l/ei,r, 1/r) of Lemma 8.1, notice that since the roots of (8.8) are single or at most dou- 
ble 18 , there are N/2 different solutions because of the reciprocity of (8.7) and (8.8). For these 
N/2 solutions (i.e., states) e\e2ese4 = —1 , and Q is given by 

Q(z) = (zy- l)(z - y) {z 2 - (t + l/t) z + l) / z 2 , (8.14) 

with 

1 2 A (r + 1/r) -4 . . 

t + - = T a > A^±l. (8.15) 

t r + l/r-2A ^ v ; 

The number A(w) is obtained inserting (8.14) and (8.15) into (5.3). This result shows also 
that (8.12) and (8.14) are more accurate trial functions to solve (5.3) than the general (8.1). 
Contrary to what we have done along this paper, we do not write the function A(v) associated 
to (8.14) and (8.15) for general A" , but we write it when A^ = 8 , which is 

A{v) = 2 a 4 6 4 + c 2 (2 A 2 a 3 b 3 - A 3 a 2 b 2 (a 2 +b 2 )-2 X l ab (a 4 + 6 4 - a 2 b 2 ) - a 6 - b 6 ) , (8.16) 

with Ai, A2, A3 certain numbers depending on A that we do not specify. The object to present 
(8.16) is to comment about the excluded cases A = ±1 pointed in (8.15). We have excluded 
these two points for mathematical reasons only. Let us fix A = 1 (we center the discussion 
in this value because the polynomial (8.8) indicates that the situation when A = —1 is the 
same just negating r ). Subtituting A = 1 in (8.15), the r.h.s. reduces either to the constant 
— 2 or to the indetermination 0/0 19 : which is then the function (8.14) and how many of them 
can one write when A = 1 ? We wont be more explicit in this point now, however we want 
to convince the reader that for A^ = 8, A = 1 there are four (eventually A^/2 for general A^ , 
if things go as they shall) bound pair states with eie2e3e4 = —1 : we have just constructed 
the states (4.1) with (4.2), (4.4) and (4.5) imposing the conditions (6.5) and (6.6); we have 
obtained exactly four states, and have checked (diagonalizing numerically the matrix block) 
that they are eigenvectos of the transfer matrix (2.1) when N = 8,n = A . The associated 
eigenvalues are precisely (8.16) with Ai = 0, -3.69963, -1.76088, 0.460505 20 , and A 2 , A 3 given 

18 The discriminant of (8.8) in r vanishes only for A = ±1/2, ±1 , thus indicating multiplicity of the roots r 
more than 1 only for these values. Why for these values? Notice that the bilinear transformation e% — » 
in the r.h.s. of (8.2) (and in the r.h.s. of (8.3) for e<t ) collapses to a constant when A = ±1/2 instead of being 
a one-to-one mapping. This justifies the multiplicities at A = ±1/2 . A similar reason happens when 6364 = 1 
and A = ±1 to the second factor in the r.h.s. of equations (8.2), (8.3) 

19 r — 1 is solution of (8.8) when A = 1 

20 Approximated to the nearest six digit number the last three data 
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in terms of Ai by 

A 2 = 2 "^r 4Al , A 3 = 2A 2 + 2A 1 -1, A = l. (8.17) 



In conclusion, for each real value of A in the vertex model, there are N '/2 + 1 bound pair 
states in the n = 4 block of the N -site transfer matrix. The number of such states is correct 21 
because exact diagonalization of the block corroborates it: our numerical experiments carried 
up to N = 12 with different but arbitrary values of the activities a, b, c confirm that the 
numbers A(v) obtained substituting Q by (8.14) with (8.15) and (8.8) into (5.3) are true 
eigenvalues of the transfer matrix. The number (8.13) is also an eigenvalue. We have no reason 
then to doubt that they are eigenvalues for general N as well. The author thus admits the 
number N/2 + 1 as absolutely right. 



For n = 5 , we count a total of N bound pairs. This is so because the solutions (e^, e±, e*>) = 
(1, 1, 1), (—1, —1, —1) of equations (8.4)-(8.6) 22 do not afford eigenvalues of the transfer matrix 
for A generic. We noticed this fact from our numerical tests carried with different values of 
a, b, c and N = 10, 12 : the numbers A obtained with (5.3) and Q as in (5.2) with zeros at 
z\ = y,Z2 = l/y,zs = Z4 = z$ = +1 and y arbitrary, do not correspond to eigenvalues of 
the transfer matrix 23 . Unlike this, the solutions in Lemma 8.2 that satisfy 636465 = — 1 afford 
N/2 bound pairs for each A , and the solutions that satisfy 636465 = 1 afford another N/2 
bound pairs (even for A = ±1 in both cases). The corresponding numbers A were checked 
numerically. These eigenvalues are obtained with 

Q(z) = (zy- l)(z - y) (z 2 - (t + 1/t) z + l) (z + l)/z b '\ (8.18) 

the plus sign in + is for 6162636465 = 1 (i.e., 636465 = —1 ), the minus sign for 6462636465 = 
— 1 . In both functions written in (8.18) 

1 2A(r + l/r)+4 . , 

t + - = V . ' J Z , A + ±1, 8.19 

t r + l/r + 2A r y ' 

but r is the root of different polynomials, as stated in Lemma 8.2. 

We write an example for = 10 and A = — 1/2 , with the choice of z, y, p as in the left 
column of Table 1. After diagonalizing numerically the blocks n = 4, 5 of the transfer matrix, 
the eigenvalues corresponding to bound pairs (we have recognized them because they match 
exactly our predicted values) approximated to the nearest six digit number are 24 : 

21 In reference [1] were found 5 states when N — 8, n — 4 , as we mentioned in Sect. 1. Our result agrees with 
that number 

22 We mentioned these solutions in the paragraph before Lemma 8.1 

23 The reason is that the states derived from these solutions proceeding as in Sections 4 and 8 are the zero 
vector 

24 An interesting question is if they can be recognized in another manner 
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n = 4 


[deg] 


n — 5 


[deg] 








0.111240 + 


[2] 


0.223 385 + 


[2] 


-0.065 746 4" 


[2] (r 


= 2) 


-0.040 142 3 " 


[2] (r = -l) 


0.138 383 + 


[2] 


-0.086 922 0" 


[2] (r 




-0.047 993 4" 


[1] 


0.074 747 8 + 


[2] 


-0.116 844" 


[1] 




-0.080 403 8" 


[1] 


0.043 830 8 + 


[2] 


-0.226 567" 


[1] 




-0.158 869" 


[1] 


0.032 535 7 + 


[2] 


-0.463 488" 


[1] 




-0.280 662 " 


[1] 













Table 2. Each eigenvalue listed is followed by a sign + or — : the sign + indicates that e^e^ = 1 (or that 
Cie2e3e4e 5 = 1 if n = 5 ), the sign — that the product of the Bethe roots is —1. The degeneration of the 
eigenvalue is [deg]. In the column corresponding to n = 4, the number 0.111240 coincides with (8.13), and the 
remaining five values agree with the theoretical A obtained inserting (8.14) and (8.15) into (5.3). The eigenvalue 
that corresponds to r = —1 , remember that in this column r is a solution of (8.8), is degenerated. This degeneration 
is not a surprise, because it is a case in which two Bethe roots coincide ( ej, = e4 = — 1 ), and when it is true that the 
eigenvector associated to such cases is usually the zero vector, when A = —1/2 it is not. Regarding the list when 
n = 5 , the values with a + correspond to solutions r of (8.10), and the values with a — to solutions r of the 
polynomial that is obtained changing in (8.10) the variables r, A by — r, — A . Totally expected is the degeneration 
of the eigenvalue —0.086 922 since e: 3 = = — l,e 5 = 1 . But the degeneration of —0.065 746 4 which happens 
for e: 3 = — 2, = — 1/2, es = —1 is less expected. 

The last comment of the paper: the numerators of (6.1), (7.5), (8.12), (8.14) and (8.18) are 
polynomials in in z with a reciprocal property, that is if R is a solution, so it is 1/R . When 
looking for other Q's in n = 7 (say) one has to restrict to numerators with this property. 
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